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Summary 

Streptomyces clavuligerus ATCC 27064 and 
S. clavuligerus AccaR::tsr cultures were grown in 
asparagine-starch medium, and samples were taken 
in the exponential and stationary growth phases. 
Transcriptomic analysis showed that the expression 
of 186 genes was altered in the ccaR-deleted mutant. 
These genes belong to the cephamycin C gene 
cluster, clavulanic acid gene cluster, clavams, 
holomycin, differentiation, carbon, nitrogen, amino 
acids or phosphate metabolism and energy produc- 
tion. All the clavulanic acid biosynthesis genes 
showed M c values in the order of -4.23. The blip gene- 
encoding a (3-lactamase inhibitory protein was also 
controlled by the cephamycin C-clavulanic acid 
cluster regulator (M c -2.54). The expression of the 
cephamycin C biosynthesis genes was greatly 
reduced in the mutant (M c values up to -7.1), while the 
genes involved in putative (3-lactam resistance were 
less affected (M c average -0.88). Genes for holomycin 
biosynthesis were upregulated. In addition, the lack 
of clavulanic acid and cephamycin production nega- 
tively affected the expression of genes for the 
clavulanic acid precursor arginine and of miscella- 
neous genes involved in nitrogen metabolism (amtB, 
glnB, glnA3, glnA2, glnA1). The transcriptomic results 
were validated by quantative reverse transcription 
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polymerase chain reaction and luciferase assay of 
/uxAB-coupled promoters. Transcriptomic analysis of 
the homologous genes of S. coelicolor validated the 
results obtained for S. clavuligerus primary metabo- 
lism genes. 

Introduction 

Streptomyces species are the largest group of antibiotic- 
producing micro-organisms. Antibiotic biosynthesis genes 
are usually clustered and their expression is regulated in 
response to nutritional environment, culture time and 
cell density. Most of the gene clusters for antibiotic 
biosynthesis contain regulators responding to these 
parameters (Martin and Liras, 2010). Regulators of LysR- 
(Perez-Redondo era/., 1998) or large ATP-binding regu- 
lators of the LuxR family (LAL)-type (Wilson era/., 2001; 
Anton era/., 2004; Tahlan etal, 2007) are present in 
several clusters, but the most common ones are the acti- 
vator proteins belonging to the SARP (Streptomyces anti- 
biotic regulatory proteins) family. They are present in 
many antibiotic clusters (Bate etal., 2002; Garg and 
Parry, 2010; He etal., 2010) and bind specific 
deoxyribonucleic acid sequences activating expression of 
antibiotic biosynthesis genes. 

Cephamycin C biosynthesis in Streptomyces cattleya is 
controlled by a SARP-type protein encoded by thnU, a 
gene that is not located in the cephamycin C cluster 
(Rodriguez etal., 2008). In Streptomyces clavuligerus, a 
similar protein, cephamycin C-clavulanic acid cluster 
regulator (CcaR), encoded by a gene located in the 
cephamycin C cluster, controls both cephamycin and 
clavulanic acid (CA) biosynthesis (Perez-Llarena etal., 
1997; Santamarta etal., 2011). The two clusters are 
located side by side in the chromosome. 

Specific sequences for the binding of SARP proteins 
have been described in antibiotic biosynthesis genes 
(Wietzorrek and Bibb, 1997). In S. clavuligerus, triple 
heptameric conserved sequences are responsible for 
CcaR binding and controls the expression of the lat, cefF, 
cefD and ccaR genes in the cephamycin C cluster. In the 
CA cluster, CcaR binds sequences upstream of ceaS2, 
which encodes the first enzyme of the CA pathway, and of 
claR, for the LysR-type regulator controlling late steps in 
CA biosynthesis (Santamarta etal., 2011). 

At the beginning of this study, the S. clavuligerus 
genome was not published, but many S. clavuligerus 
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genes had been deposited in the database and addi- 
tional genes were available through the DSM (Delft, The 
Netherlands) S. clavuligerus sequencing project. There- 
fore, we constructed microarrays containing probes for 
800 genes of the S. clavuligerus genome as well as for 
7728 genes of the Streptomyces coelicolor genome in 
order to compare the housekeeping genes of both strains. 
This microarray has been used to detect gene expression 
in different S. clavuligerus mutants. In this article, we 
report the results observed for S. clavuligerus AccaRr.tsr, 
a strain lacking the SARP regulator that controls 
cephamycin C and CA biosynthesis. 

Results 

Transcriptomic analysis of antibiotic biosynthesis 
genes in S. clavuligerus ATCC 27064 and 
S. clavuligerus 4ccaR::tsr 

Streptomyces clavuligerus ATCC 27064 and 
S. clavuligerus AccaR::tsr cultures were grown in 
asparagine-starch (SA) medium, and samples were taken 
in the exponential and stationary growth phases (Fig. 1) 
as indicated in Experimental Procedures. The microarrays 
analysis showed differences in expression in 186 genes 
of the ccaff-negative mutant at one specific sampling time 
at least (Table 1). 

These genes belong to the cephamycin C gene cluster, 
CA gene cluster, clavams, holomycin, cellular differentia- 
tion, carbon, nitrogen, amino acids or phosphate metabo- 
lism and energy production. 

CA biosynthesis genes. Transcriptomic studies showed 
that expression of the four genes of the ceaS2 to cas2 
operon, encoding enzymes for the early steps of CA 
biosynthesis, dropped in the mutant to an average M c 
value of -6.34 in the exponential phase. Expression of 
claR, encoding the regulatory protein ClaR, is under CcaR 
control (Perez-Redondo era/., 1998; Santamarta era/., 
2011). This gene was underexpressed in the ccaR- 
deleted mutant (M c -4.61 ) and concomitantly transcription 
of the genes for the late steps of the pathway car, gcaS as 
well as cyp, orf12, orf13, orf14, oppA2, oppA1 and orf16, 
of unknown function, dropped an average M c value of 
-4.23, being oat2, orfl 3 and car above this value. 

Expression of orf18 and orf19, which encode penicillin- 
binding proteins, orf20 for a cytochrome P450 and pcbR, 
tentatively involved in p-lactam resistance, was lower in 
the mutant than in the wild-type strain. All these genes 
were expressed at higher M c values (average M c -0.88) 
than the biosynthetic genes in the ccafl-mutant. 

An interesting finding was made in relation to genes 
encoding (3-lactamase proteins. The blp gene, encoding a 
putative |3-lactamase inhibitory protein, which is located 
downstream of ccaR, was clearly underexpressed in the 
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Fig. 1. Growth and antibiotic production in Asparagine-starch 
medium. Strains used: S. clavuligerus ATCC 27064 (open symbols) 
and S. clavuligerus AccaRr.tsr (closed symbols). A. Growth (circles) 
as measured by the DNA content. B. Cephamycin C (circles) or 
clavulanic acid (triangles) production. C. Holomycin production 
(squares). 

mutant (M c -6.2 and -4.8 in the exponential and station- 
ary phases), which might be due to the absence of the 
ccaR promoter in the ccafi-mutant. The second gene, 
blip, encoding a well-characterized (3-lactamase inhibitory 
protein (Doran era/., 1990), was also underexpressed in 
the mutant. This gene, located outside the CA gene 
cluster, presented M c values of -2.54 and -2.19 in the 
exponential and stationary phases, respectively. 
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Table 1. Streptomyces clavuligerus AccaR::tsr gene expression as compared with S. clavuligerus ATCC 27064. 



Exponential phase Stationary phase 



Code 


Gene 


Product 


Mc 


FDR 


Mc 


FDR 
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Holomycin biosynthesis 
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Table 1 . cont. 



Exponential phase Stationary phase 



Code 
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Product 
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< 1 E -06 


-0.95 


/ OOI — —04 

4.30b 


bOLAv_oyby 




Fumarate reductase iron-sulfur subunit 


-1 .67 


1 .82 E 


—1.15 


1 .bbb 


O /""■» 1 All A~!C^~7 

bOLAV_4767 




Monophosphatase 


-1 .86 


, j i— -06 

< 1 b 


-0.34 


2.1 6E" 01 


Carbon metabolism 














OPI AW f\COH 

bOLAV_Ubo 1 


glpF2 


Putative glycerol uptake facilitator protein 


—1 .49 


1 .1 2E 


-1 .95 


< 1 b 


SCLAV_0632 


glpK2 


Putative glycerol kinase 


-1.13 


r\ r~r\ i — -03 

9.59E 


-1 .82 


< 1 E" 06 


bGLAV_087b 


gylR 


Glycerol operon regulatory protein 


-0.31 


O fr~il — -01 

3.62b 


-1 .34 


< 1 E" 06 


c*/~>i aw r\o~7~7 

bGLAV_0877 


glpF1 


Putative glycerol uptake facilitator protein 


0.30 


6.1 9E" 01 


-0.42 


2.76E" 01 


o /*"* i aw no70 
bOLAV_Uo/o 


gipKi 


Glycerol kinase 


-0.07 


y.Uob 


-1.19 


i .yyb 


O /""■» 1 aii no-7fi 

bOLAV_0879 


glpD 


Glycerol-3-phosphate dehydrogenase 


-0.02 


9.35E" 01 


-0.94 


< 1 E" 06 


bGLAV_5509 


gap2 


Glyceraldehyde-3-phosphate dehydrogenase 2 


-2.61 


< 1 E -06 


-1 .75 


< 1 E" 06 


o^i aw /icon 
bULAV_4o2y 


glcP 


Glucose permease 


1 .79 


< 1 b 


-0.60 


"7 (T7C-02 


O /~> 1 All «Art"7C 

bGLAV_p0975 




Ribulose-phosphate 3-epimerase 


-1 .25 


, j i— -06 

< 1 b 


-1 .36 


< 1 E" 06 


Nitrogen metabolism 














SCLAV_4534 


amtB 


Ammonium transporter 


-4.10 


< 1 E 


-5.04 


< 1 E 


O /~> i aii ^nroc: 

bOLAV_4535 


ginB 


Putative nitrogen regulatory protein P-ll 


-4.39 


< 1 E~ 06 


-5.09 


< 1 E" 06 


bULAV_p1 4b2 


giniii 


Glutamine synthetase III 


-0.60 


O.oob 


—1 .54 


< lb 


bOLAV_1431 


gin A3 


Glutamine synthetase (Glutamate-ammonia ligase) 


-4.67 


< 1 b 


-5.32 


< 1 E" 06 


OPI A\ / H/HC 


glnA2 


Glutamine synthetase I (Glutamate-ammonia ligase I) 


-3.56 


< 1 b 


-2.95 


< 1 b 


SCLAV_1473 




Glutamine synthetase 


-0.88 


6.84E" 03 


-2.35 


<1E -Q6 


C/^l A\/ ACCZC\ 

oL/LAV_4oOU 


giuui 


Glutamate transporter permease 


— U.o4 


1 .b^ib 




4.1 4b 


SCLAV_0834 


glnA1 


Putative glutamine synthetase 


-1.77 


< 1 E" 06 


-2.41 


< 1 E" 06 


Phosphate metabolism 














ooLAV_l / 1 y 


phoH 


Phosphate starvation-induced protein 


1 .25 


04 

1 .37E 


-0.06 


8.67E 


h>GLAV_31 66 


pstB 


PstB protein — Phosphate import ATP-binding 


-0.93 


A A A 1— -02 

1 .41 b 


-1 .50 


1 .14E" 05 


oULAV_o1b/ 


pstA 


PstA protein — Permease component 


—1 .12 


/I O A C-05 

4.o4b 


-1 .61 


< 1 b 


bOLAV_316o 


pstC 


PstC protein — Permease component 


-1 .11 


2.41 E" 02 


—1 .77 


4.12E" 05 


oULAV_o1 by 


pstS 


PstS protein precursor — Periplasmic component 


-0.93 


3.1 OE" 02 


—1 .41 


1 .bbb 


O^l AW OOOA 


phoU 


Putative phosphate transport system regulatory 
pruicin 


-1 .01 


O COL-04 

o.o2b 


-1 .23 


< 1 b 


Lipid metabolism 














SCLAV_4986 




Putative acetyl-coenzyme A synthetase 


-0.55 


4.01 E" 02 


-2.18 


< 1 E- 06 


SCLAV_3406 




Putative acetyl/propionyl CoA carboxylase alpha 


-1.49 


2.87E" 03 


-2.15 


< 1 E -06 


SCLAV_3405 




Putative acetyl/propionyl CoA carboxylase beta 


-1.64 


< 1 E" 06 


-2.41 


< 1 E" 06 


Transcriptional and 














regulatory proteins 














SCLAV_p0826 




Putative AraC-family transcriptional regulator 


0.88 


2.87E" 05 


-0.05 


8.17E" 01 


SCLAV_p0894 


brp 


Gamma-butyrolactone receptor protein 


0.86 


7.97E" 04 


-0.55 


1.98E" 02 


SCLAV_p1319 




Putative transcriptional regulator AraC family 


-1.72 


1 .75E" 04 


-1.10 


1.03E" 02 
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Exponential phase Stationary phase 



Code 


Gene 


Product 


Mc 


FDR 


Mc 


FDR 


o /■-« i a\ / h aaa 

oOLAV_1096 




Transcriptional regulator, GntR-family protein 


-0.66 


H A A 1 — -OA 

1 .20b 


-0.09 


6.1 6E 


f\ f \ 1 AW -I A AA 

SCLAV_1433 




Putative regulatory protein 


-1.18 


A AO 1 — —OP 

2.33E 


-0.60 


1 .94E" 01 


oULAV_l b2 1 




Putative MerR-family transcriptional regulator 


0.67 


2.0OC 


-0.93 


2.y2b 


f \ 1 Al 1 HAf 7 

SCLAV_1957 


adpA 


AraC-family transcriptional regulator 


-1 .22 


5.42E" 04 


-1 .43 


1 .52E" 05 


SCLAV_1958 


ornA 


Oligoribonuclease 


-1 .00 


< 1 E" 06 


-0.97 


< 1 E" 06 


SCLAV_2732 




Two-component transcriptional regulator 


1 .31 


2.20E" 03 


-0.23 


6.00E" 01 


f\ f \ 1 ai i A A A H 

SCLAV_3001 




Putative gntR-family transcriptional regulator 


-0.68 


A "HOI — —04- 

8.78E 


-1.51 


< 1 E~ 06 


SCLAV_4054 




WhiB-family transcriptional regulator 


2.62 


< 1 E" 06 


0.32 


3.1 3E" 01 


f\ /—^ 1 A \ 1 j1AA~7 

SCLAV_4937 




Putative regulatory protein 


1 .05 


A ~~7AI — —OR 

9.79E "° 


-1 .06 


< 1 E" 06 


O Z" 1 1 A\ / ro7n 

oGLAV_527o 




AmphRI-like transcriptional regulator 


2.01 


A tz A 1^-04 

2.59b 


1 .36 


8.01 b 


Unknown function 














Of*\ aw aah o 

oOLAV_001 o 




Cytochrome P450 monooxygenase 


-1.16 


o.47b 


—1 .47 


1 .77E 


r^z"**i aw a a a a 

SCLAV_0633 




ATP— GTP-binding protein 


-0.78 


A A A 1 — —OP 

4.44E 


-0.52 


1 . 1 9E" 01 


C^f^-l AA / AAAA 

SCLAV_0636 




Putative large secreted protein 


-0.76 


1 .59E" 02 


-0.17 


5.91 E~ 01 


SCLAV_0646 




Putative inhibitor of KinA 


-0.93 


A A A 1 — —02 

3.99E 


0.08 


8.53E" 01 


O Z" 1 1 AW A~7 ,1 A 

oGLAV_0743 




Peroxidase 


1 .46 


< 1 E" 06 


1 .40 


< 1 E" 06 


r^z"**i aw h a a r~ 

SCLAV_1335 




Two-component system sensor kinase 


0.51 


3.31 E~ 02 


-0.37 


7.66E~ 02 


O Z" 1 1 A\ / -\ >1 A A 

oOLAV_1344 




Conserved phosphoesterase 


-0.78 


1 .36E~ 02 


—1 .11 


O A H I^-OR 

8.41 b 


CPI AW H CC/I 

oOLAv_ 1 Oo4 




Acetyl-CoA acetyltransferase 


-2.06 


< 1 E" 06 


-2.68 


,. h r--06 

< l b 


r^z"**l AW n r^r 

SCLAV_1565 




Cytochrome P450 hydroxylase 


-1 .26 


< 1 E~ 06 


-2.67 


< 1 E~ 06 


CPI AW 




Hypothetical protein 


1 .02 


<lb 


0.03 


8.67E -01 


f\ /"""•» 1 AW H ~7 A 

SCLAV_1748 




DUF143 domain-containing protein 


0.45 


1 .28E" 02 


0.14 


4.26E" 01 


f\ Z~^ 1 A \ 1 H A f - A 

SCLAV_1959 




Sensor protein 


-0.91 


< 1 E"°6 


-1 .36 


< 1 E -06 


f\ 1 AW AAAA 

SCLAV_2623 




SclavP3 predicted orf 


1 .03 


1 .10E~ 02 


0.05 


9.07E -01 


r^^i aa / a a a r~ 

SCLAV_2625 




Subtilase-type protease inhibitor precursor 


0.96 


A A A 1 — —OP 

4.44b ^ 


-0.04 


A AAI — —01 

9.33E U1 


C/" 1 ! A\ / O i ("1/1 

oULAV_o1 y4 




DUF1416 domain-containing protein 


-0.77 


2.64E" 02 


0.55 


7.31 E" 02 


O Z" 1 1 A\ / A H AH 

oOLAV_413 1 




Metallophosphoesterase 


-0.80 


3.06b 


-1 .40 


1 .54E" 05 


f \ 1 AA / /1AAA 

SCLAV_4308 




Methylmalonyl-CoA epimerase 


-0.68 


1 .77E" 05 


-1 .32 


1 .31 E~ 15 


O z™» 1 A\ / /I A C A 

oOLAV_4352 




Integrin-like protein 


1.15 


7.63E" 03 


0.42 


A AAI — -01 

2.9ob 


f\ s~\ i aa / iiorr 

SCLAV_4355 




Hypothetical protein 


-0.64 


A AAI — —OP 

2.99b Uii 


-1 .35 


< 1 E" 06 


cpi aw /iocn 
oOLAV_4oo9 




Metalloendo peptidase 


-2.95 


1 .92b 


-3.36 


, h r--06 

< i b 


oCLAV_4530 




Acetiltransferase 


2.52 


< 1 b 


-0.21 


f A A 1 — -01 

6.30b 


r^/^l AA / A~7 -i ~7 

SCLAV_471 7 




Putative hydroxylase 


-1 .20 


A AAI — —01 

6.98E UJ 


-1 .71 


2.05E" 05 


O Z" 1 1 A\ / C~\A<~\ 

oOLAV_5249 




Membrane protein 


1 .31 


A A A 1^-04 

6.40b 


0.20 


6.1 0E~ 01 


oOLAV_p0763 




Amidohydrolase:Amidohydrolase-like precursor 


-0.84 


3.63b ^ 


-0.1 0 


7.94E" 01 


f\ 1 AA / -~ H H A A 

SCLAV_p1123 




Putative methyltransferase 


3.26 


< 1 E" 06 


3.44 


< 1 E" 06 


oOLAV_p1142 




Rhs family protein 


-1 .01 


< 1 E" 06 


—1 .77 


< 1 E" 06 


o i a nun™ :—i 

pbuL2 rlasmid 














C*«lrtAO AH AA A A~7 A Ar 

bClaA2_0101 00027605 




Helicase 


1 .02 


O A A [Z-03 

3.06b 


-0.11 


01 

7.59E 


r~\ i a r*i ahahaaaa - 7 a h a 

SclaA2_01 01 00027610 




Hypothetical protein 


0.91 


a a r~ 1 — —04 

3.35E 


0.14 


6.04E" 01 


C*«lmAO ah ah AAnnroc: 

oClaA2_0101 00027625 




Hypothetical protein 


0.91 


A A A 1 — -OR 

4.42b 


1 .27 


< 1 b 


C*«lrtAO A H A H AAAA~7AAA 

bClaA2_0101 00027920 




Transposase 


1 .00 


7.48E" 04 


0.11 


7.1 6E~ 01 


[~\ 1 a A AH AH AAAA~7AAA 

SclaA2_0101 00027930 




Hypothetical protein 


2.15 


< 1 E" 06 


1 .27 


1 .21 E" 05 


[~\ 1 a A AH AH AAAA - 7AAI~~ 

SclaA2_0101 00027935 




Hypothetical protein 


1 .07 


A A A 1 — —03 

4.80E 


0.86 


1 .27E~ 02 


C*«lrtAA AHAH A A A A~7 AIT ET 

bClaA2_0101 00027955 




GntR-family regulatory protein 


0.85 


A OA 1 — -09 

3.89fc 


0.69 


5.42E" 02 


OrtlrtAA AHAH A A A A~7 A~7ET 

oClaA2_0101 00027975 




Hypothetical protein 


0.50 


2.02b 


0.96 


^ j 1— -06 

< 1 b 


f-\ 1 A A AH AH AAAA - 7AAA 

SclaA2_0101 00027990 




Hypothetical protein 


0.70 


6.76E" 03 


0.11 


6.74E" 01 


oClaA2_0101 00028020 




Hypothetical protein 


1 .64 


8.91 b 


1.18 


5.1 5E* 05 


oClaA2_01 01 0002801 5 




Phosphatase 


0.88 


8.96b 


-0.1 6 


f AAI — -01 

6.20b 


OClaA2_U1 01 0002oo2o 




Hypothetical protein 


1 .85 


< 1 b 


0.42 


1 .64E" 01 


oClaA2_0101 00028330 




Hypothetical protein 


2.71 


, j 1— -06 

< 1 b 


0.98 


4.07E~ 02 


SclaA2_01 01 00028335 




Transferase 


2.97 


< 1 E -06 


1.82 


3.1 9E" 05 


SclaA2_01 01 00028340 




Telomere-associated protein 


1.15 


6.58E" 03 


0.27 


5.1 2E" 01 


SclaA2_01 01 00028350 




Hypothetical protein 


0.73 


6.58E" 03 


-0.29 


2.42E" 01 


SclaA2_01 01 00028360 




Hypothetical protein 


0.75 


3.28E" 02 


0.10 


7.66E" 01 


SclaA2_01 01 00028340 




Telomere-associated protein 


0.79 


1.61E" 02 


0.02 


9.56E" 01 


pSCL1 Plasmid 














SclaA2_01 01 00027570 




Hypothetical protein 


1.69 


4.80E- 03 


0.24 


6.85E- 01 


SclaA2_01 01 00027560 




Hypothetical protein 


1.95 


4.33E" 13 


1.07 


1 .82E" 05 


SclaA2_01 01 00027550 




Hypothetical protein 


1.11 


1.81E" 02 


-0.22 


6.21 E" 01 


SclaA2_01 01 00027545 




Hypothetical protein 


1.11 


6.58E- 03 


-0.51 


1 .73E" 01 



ABC, ATP-binding cassette transporters; ATP, adenosine triphosphate; BLIP, p-lactamase-inhibitory protein; GTP, guanosine triphosphate; NADH, 
nicotinamide adenine dinucleotide. 
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A gene external to the CA cluster, adpA, encoding a 
regulatory protein involved in antibiotic biosynthesis 
(Lopez-Garcia era/., 2010) was also slightly 
downregulated. 

Genes for CA precursors. Clavulanic acid biosyn- 
thesis occurs through condensation of arginine and 
glyceraldehyde-3-phosphate by the carboxyethylarginine- 
synthase, the first enzyme of the CA pathway. 
Streptomyces clavuligerus and S. coelicolor probes of 
genes for arginine biosynthesis and glycerol utilization 
were present in the microarrays. 

The arginine biosynthesis genes were downregulated 
in the ccaff-deleted mutant. Six genes for arginine 
biosynthesis, located in three separated locations in the 
S. clavuligerus genome, showed average M c values of 
-2.01 and -2.85 in the exponential and stationary phases 
respectively. These genes are controlled by the argR- 
encoded regulator ArgR, whose expression decreased to 
-1.45 and -2.34, at both sampling times. 

The co-transcribed genes for glycerol utilization, glpF1- 
glpK1-glpD1, were weakly but significantly affected by the 
lack of the CcaR regulator, and the same occurred to the 
gylR gene, encoding the glycerol utilization regulatory 
protein (Bahos era/., 2009). 

However, a positive effect of CcaR on the second 
cluster for glycerol utilization might occur, because glpF2 
and glpK2 were underexpressed in the ccaff-deleted 
strain. The most affected gene in this group was gap2, 
encoding glyceraldehyde-3-phosphate dehydrogenase 
(Mc -2.61 in the exponential growth phase). 

Cephamycin C biosynthesis genes. In relation to 
cephamycin C biosynthesis, the transcriptomic data fitted 
well with those previously obtained by quantative reverse 
transcription polymerase chain reaction (Santamarta 
era/., 2011). The M c values for lat, cmcl, cmcJ, cefD and 
cefF genes, either starting transcriptional units or close to 
the first gene of the operon, showed values ranging from 
-6.15 to -7.48 in the exponential phase, close to the 
-7.12 value of ccaR (which is deleted in the mutant), 
indicating a complete lack of transcription of these units. 
All other genes related to cephamycin biosynthesis 
showed Mc values between -3.2 and -7.12 in the expo- 
nential phase and their expression dropped in the ccaff- 
mutant by an average of 18% in the stationary phase. 

The pbp74 and bla genes were barely affected in the 
exponential phase but were less expressed in the mutant 
in the stationary phase. 

Genes for other antibiotics produced by S. clavuligerus. 
All the holomycin biosynthesis genes tested were 
overexpressed in the ccaff-deleted mutant, with an 
average M c value of 5.57 in the exponential phase, while 



the regulatory gene hlmM had a M c value of 2.01 . These 
results coincide with the formation of holomycin by the 
S. clavuligerus ccaR::aph disrupted mutant detected by 
Fuente and colleagues (2002). A downregulation was 
observed in the clavams and CA paralogous clusters, 
especially in the stationary phase. Production of clavams is 
relatively variable and medium-dependent. However, cas1 
encoding the second clavaminate synthase isoenzyme 
(M c -3.32 and -4.29) and cvm8, encoding a serine 
hydroxymethyltransferase (M c -2.59 and -3.18) were 
strongly downregulated in the ccaff-deleted mutant. 

Nitrogen and phosphate metabolism 

The lack of CcaR and, therefore, of CA and cephamycin C 
formation, affected nitrogen metabolism. A strong 
downregulation was observed in the exponential phase in 
the expression of the ammonium transporter-encoding 
gene amtB (M c -4.1 ) and the downstream gene, glnB (M c 
-4.39), which encodes the uridilyltransferase regulatory 
protein Pll. This effect was stronger in the stationary 
phase (M c -5.04 and -5.09, for amtB and glnB respec- 
tively). Also genes for glutamine synthetases {glnA3, 
glnA2, glnA1) were underexpressed in the ccaff-deleted 
mutant (M c -4.67, -3.56 and -1.77) in the exponential 
phase. A weaker, but clear downregulation also occurred 
in g/uD1, encoding a glutamate permease. The homolo- 
gous probes of S. coelicolor (also present in the arrays) 
for glutamine synthetase {glnA, glnll, SC01 61 3) as well of 
genes for urease (ureAB) and an allantoinase (SC06247) 
gave a lower signal when hybridized with S. clavuligerus 
AccaR::tsr messenger ribonucleic acid (mRNA) in relation 
to S. clavuligerus ATCC 27064 mRNA. This might indicate 
that cephamycin C and CA production requires a strong 
demand of nitrogen-derived precursors and, therefore, 
nitrogen metabolism slows down when antibiotic produc- 
tion is blocked. 

Some genes involved in phosphate transport {pstA, 
pstCandpstS) showed a weak but clear underexpression 
in the ccaff-mutant, especially in the stationary phase 
sampling point. 

Miscellaneous genes 

Several genes involved in energy production or lipid 
metabolism were affected in the ccaff-deleted strain. This 
was the case for the genes encoding putative acetyl-CoA 
synthetases and acetyl-CoA carboxylases (SCLAV_4986, 
SCLAV_3405, SCLAV_3406) and the same occurred with 
the acetyltransferase- and methyltransferase-encoding 
genes (SCLAV_4530 and SCLAV_p1123). A set of genes 
whose function is unknown (SclaA2_01 01 00027930, 
SclaA2_01 01 00028330 and SclaA2_01 01 00028335) all 
located in plasmid pSCL2, showed a clear overexpression 
in the mutant in the exponential phase of growth. 
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Effect on regulatory genes. The transcription of 1 3 genes 
encoding different types of regulators was affected in the 
mutant, including adpA (see previous discussion). A very 
strong downregulation was observed for the wft/'B-family 
transcriptional regulator (M c 2.62) and for the gene encod- 
ing an AmphR1-like regulator (SCLAV_5278) with an M c 
value of 2.0 in the exponential phase. 

Expression of S. clavuligerus genes as detected 
by transcriptomic analysis of orthologous 
S. coelicolor genes 

The whole S. coelicolor genome was represented in the 
microarrays. Therefore, we were able to compare expres- 
sion of miscellaneous genes by heterologous hybridiza- 
tion of S. coelicolor probes present in the microarray with 
mRNA from S. clavuligerus. Heterologous transcriptomic 
studies confirmed the lower expression of the glnA, glnll, 
SCLAV_0834, ureAB and gap2 genes in the AccaR 
mutant. Other blocks of underexpressed genes were 
those encoding succinate dehydrogenases (SCLAV_3969 
and SCLAV_3970) and the nuo genes (nuoEFJKLMN) 
for different subunits of the nicotinamide adenine 
dinucleotide (NADH) dehydrogenase. Genes for mem- 
brane proteins or for hypothetical proteins were also 
found to be underexpressed; but some were over- 
expressed, including the orthologous aminotransferase 
encoded by SCLAV_5663 (M c 2.63) and the 
methyltransferase encoded by SCLAV_5654 (M c 2.56) 
(Supporting Information Table S1). 

Validation studies 

Quantitative reverse transcription polymerase chain 
reaction (qRT-PCR). The validation of many genes differ- 
entially expressed in this transcriptomic analysis has 
already been carried out using proteomics or qRT-PCR 
(Santamarta etal., 2011). Such is the case for all the 
downregulated genes of the CA pathway and for the lat, 
cmcl, cefD or cmcT genes in the cephamycin C gene 
cluster. 

Other differentially regulated genes present in Table 1 
and Supporting Information Table S1 were validated 
by qRT-PCR performed on reverse-transcribed RNA 
samples. A total of 15 genes (pcbC, pcd, claR, ceaS2, 
blip, html, hlmA, gapA2, glnA1, glnA2, SCLAV_4359, 
SCLAV_3668, mprA2, SCLAV_5661 and, as control, hrdB) 
were chosen for analysis (Fig. 2). The Pearson's correla- 
tion coefficient (R 2 = 0.9243) between M c values and rela- 
tive expression values of qRT-PCR indicates a good 
validation of the results. 

Luciferase activity. Promoters of the argR, phoH, amtB 
and SCLAV_p1123 genes were coupled to the luxAB 
genes and the constructions were introduced in S. 



clavuligerus ATCC 27064 and S. clavuligerus AccaR::tsr. 
Luciferase activity was lower in S. clavuligerus AccaR::tsr 
when the luxAB genes were expressed from the argR and 
phoH promoters (Fig. 2C) and higher when expression 
was carried out from the amtB and SCLAV_p1123 pro- 
moters in agreement with the results obtained in 
the transcriptomic experiments. The correlation index 
between the data obtained, measuring luciferase activity 
and the M c values of the respective genes, was 0.796. 
Luciferase activity differed slightly from the expected data 
of the microarrays when the enzyme was formed in the 
stationary phase using the phoH promoter. 

Discussion 

A transcriptomic study reveals the complete gene expres- 
sion in a micro-organism, allowing the comparison of over- 
producing strains and wild-type strains (Medema etal., 
201 0). However, while this method is convenient for under- 
standing the metabolic flow in industrial strains, these 
strains have been obtained through heavy mutagenesis 
and selection programs; therefore, the final results 
observed might not be targeted at a particular mutation but 
derive from indirect effects of randomly mutated genes. 

In this study, we compared the wild-type strain 
S. clavuligerus ATCC 27064 and a ccaf?-deleted mutant 
obtained by molecular methods (Alexander and Jensen, 
1998). As expected, the CA and cephamycin C 
biosynthesis genes are strongly downregulated in the 
mutant in agreement with the activator effect found previ- 
ously for CcaR (Perez-Llarena etal., 1997; Santamarta 
etal., 2011). These differences are very marked (a 37-fold 
decrease in the exponential phase) for all the CA 
biosynthesis genes, or the cephamycin biosynthesis 
genes (559-fold) but lower (average 1 .6 fold) for genes 
located at the ends of the clusters (orf18, orfW, orf20, 
pcbR, pbpA), which are supposed to be involved in 
(3-lactam antibiotic resistance. Interestingly, two genes 
encoding p-lactamase inhibitory proteins, blip and blp, are 
downregulated in the mutant. This is the first description of 
a relationship between CcaR and blip expression, a 
gene located outside the CA cluster. According to the 
transcriptomic data, blip is expressed as an independent 
gene. A putative heptameric sequence for CcaR binding is 
present 104 bp upstream of blip, suggesting a direct effect 
of CcaR on the expression of this gene; however, electro- 
phoretic mobility shift assay analysis has to be performed 
to confirm this hypothesis. The function of the blp- 
encoded protein has not been clearly elucidated (Thai 
etal., 2001), but this gene is strongly downregulated in 
the ccafl-mutant, perhaps due to a co-transcription with 
the ccaR gene located upstream. 

Medema and colleagues (2010) found several 
uncharacterized gene clusters, including one for 
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Fig. 2. Validation of the microarray results by qRT-PCR. Quantitative RT-PCR of genes indicated below using the oligonucleotides shown 
in Supporting Information Table S2. The relative values are referred to 1, assigned relative value for the expression of each gene in 
S. clavuligerus ATCC 27064. Error bars were calculated by measuring the standard deviation among biological replicates of each sample. 
A. Antibiotic biosynthesis genes. B. Miscellaneous genes indicated below. C. Validation of the microarray data using the luciferase coupled 
method. Streptomyces clavuligerus ATCC 27064 (black bars); S. clavuligerus AccaR::tsr (white bars). The numbers indicate for each gene the 
relative M c values of S. clavuligerus AccaR::tsr for each gene. 
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secondary metabolism, that were overexpressed in 
the high-producing strain. The opposite is found in S. 
clavuligerus AccaRv.tsr, in which all the genes for 
holomycin biosynthesis (hlmA to hlmM) are over- 
expressed (average values of 32-fold) confirming that 
the ccafi-knockout mutant produces more holomycin than 
the wild-type strain (Fuente era/., 2002; Robles-Reglero 
era/., 2013). These results agree with those of Chen and 
colleagues (2009) who reported an overexpression of the 
filipin biosynthesis genes in an ave7-disrupted mutant 
unable to produce avermectin and suggest that precursor 
and energy flow might be directed to the production of a 
different secondary metabolite once the most abundant 
metabolite pathway has been disrupted. 

Expression of the clavams biosynthesis cluster or the 
'paralogous' CA cluster, including the cvm7 gene, encod- 
ing a LAL regulator, was weakly downregulated except for 
the cas1 and orf8 genes, encoding the clavaminate 
synthetasel isoenzyme and a serine hydroxymethyl- 
transferase, which were strongly downregulated. This 
result contrasts with that of Medema and colleagues 
(2010) who described a significant overexpression of all 
these genes in the CA overproducing mutant that 
overexpresses ccaR. 

The effects of CcaR absence on the transcription of 
S. clavuligerus primary metabolism genes are validated 
by the parallel transcription results for S. coelicolor genes 
present in the microarray. 

Clavulanic acid is formed from arginine and a C3 com- 
pound derived from glycerol. All genes for arginine 
biosynthesis are strongly downregulated in the mutant 
suggesting that cells detect the lack of arginine require- 
ment for CA biosynthesis in the ccafl-deleted mutant and 
shut down the arginine pathway. Genes for glycerol utili- 
zation (glpF2, glpK2) are expressed in the CA non- 
producing mutants in the order of 0.4-fold in relation to the 
control strain, which is a small reduction if compared with 
the two-fold expression increase observed for these 
genes by Medema and colleagues (2010) in the high 
CA-producing strain; however, our result in the gap2 gene, 
for glyceraldehyde-3-phosphate dehydrogenase 2, indi- 
cates that this gene is downregulated (0.16-fold), confirm- 
ing what was observed by Medema and colleagues (201 0) 
in the high CA-producing strain and the relevance of gap2- 
disruption in increasing the glyceraldehyde-3-phoshate 
required for CA production (Li and Townsend, 2006). 

In summary, our results confirm most of those previ- 
ously obtained for a CcaR overproducing strain. CcaR 
binds heptameric sequences in many CAand cephamycin 
C genes. However, no clear heptameric sequences have 
been found in nitrogen metabolism genes, genes control- 
ling energy flow or genes for the antibiotic precursors 
suggesting that the lack of CA and cephamycin C produc- 
tion directly affects the flow of these pathways. 



Experimental procedures 

Culture conditions 

Streptomyces clavuligerus ATCC 27064 and the mutants 
S. clavuligerus AccaR::tsr (Alexander and Jensen, 1998) and 
S. clavuligerus ccaRr.aph (Perez-Llarena era/., 1997) were 
used in this work. Strain S. clavuligerus AccaRy.tsr was 
chosen for the transcriptomic experiments because the ccaR 
was deleted. To determine the more homogeneous and 
repetitive conditions for RNA sampling, the wild-type strain 
and S. clavuligerus AccaR::tsrwere grown in SA medium and 
DNA, cell dry weight and antibiotics production were 
analyzed. Trypticase soy broth medium (100 ml) was inocu- 
lated with 1 ml of frozen mycelia, and the culture was grown 
to an OD 600 nm of 6.5. This culture was used to inoculate (5% 
v/v) 500 ml baffled flasks containing 100 ml of semidefined 
SA medium (Aidoo etal., 1994). The cultures were main- 
tained at 28°C with 220 r.p.m. shaking. Exponential phase 
sampling was done when the cultures reached a DNA content 
of 75 and 80 |xg mT 1 in the wild type and the mutant (24 and 
32 h respectively) and 160 and 170 |xg DNA mh 1 in the wild 
type and the mutant for the early stationary phase (40 and 
50 h respectively) (Fig. 1). 



RNA isolation and purification 

Samples (2 ml) from S. clavuligerus ATCC 27064 and 
S. clavuligerus AccaR::tsr in the exponential and stationary 
growth phase in SA medium were stabilized with two volumes 
of RNA Protect Bacteria Reagent (Qiagen) for 5 min, then 1% 
(3-mercaptoethanol was added. The samples were treated as 
indicated by Alvarez-Alvarez and colleagues (2013). 



Labelling and microarray hybridizations 

Streptomyces clavuligerus microarrays were obtained from 
Agilent Technologies (Santa Clara, CA, USA), in the Agilent 
8x15K format. They include S. clavuligerus quadruple 
probes for about 800 genes and intergenic regions of some 
clusters involved in secondary metabolism, and also dupli- 
cated probes for 7728 chromosomal genes (out of 7825) of 
the S. coelicolor genome. Four biological replicates were 
made for each condition (two strains and two growth times). 
Labelling reactions were performed according to the recom- 
mendations described by BioPrime Array CGH Genomic 
Labelling Systems (Life Technologies, Carlsbad, CA, USA). 
Total RNA was labelled with Cy3-dCTP (Amersham, Freiburg, 
Germany) using random primers and Superscript II reverse 
transcriptase (Invitrogen, Carlsbad, CA, USA). Genomic 
deoxyribonucleic acid (gDNA) was labelled with Cy5-dCTP 
(Amersham) from random primers extended with the Klenow 
fragment of DNA polymerase (Roche, Basel, Switzerland). 
The final products were purified with MinElute columns 
(Qiagen, Venlo, the Netherlands) and labelling efficiencies 
were quantified spectrophotometrically. Cy3-cDNA (300 ng) 
and Cy5-labelled gDNA (10 pmol) were mixed, vacuum dried, 
resuspended in 32 p:L of hybridization solution (Agilent) and 
applied on the microarray surface. Hybridizations were 
carried out at 55°C and extended to 60 h to improve the 
quality of the results (Sartor era/., 2004). Washing, scanning 
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with an Agilent DNA Microarray Scanner G2565BA and 
image quantification were carried out as indicated previously 
(Rodriguez-Garcia etal., 2007). 

Identification of differentially transcribed genes and 
transcription profile classification 

Microarray data were normalized and analyzed with the 
Bioconductor package limma (Smyth and Speed, 2003; 
Smyth, 2004). Weighted median was applied within arrays. 
Weights were assigned as follows: 1 , probes corresponding 
to S. coelicolor genes showing a raw Cy3 intensity value 
higher than 2000; 0.25, probes corresponding to S. 
clavuligerus genes; and 0, the remaining probes. The nor- 
malized log 2 of Cy3/Cy5 intensities is referred to as the M„ 
value, which is proportional to the abundance of transcripts 
for a particular gene (Mehra etal., 2006). The procedure by 
Smyth and colleagues (2005) to include the information from 
within-array replicates was applied to the set of quadruple 
probes. The M„ transcription values of the four experimental 
conditions were compared using two contrasts, mutant 
versus wild type, corresponding to the two studied growth 
phases (exponential and stationary). For each gene, the M c 
value is the binary log of the differential transcription between 
the mutant and the wild strain. A positive M c value indicates 
upregulation, and a negative one, downregulation. False dis- 
covery rate (FDR) correction for multiple testing was applied. 
A result was considered as statistically significant if the 
FDR-corrected P-value < 0.05. The microarray data have 
been deposited in National Center for Biotechnology 
Information-Gene Expression Omnibus under accession 
number GSE51435. 

qRT-PCR 

The oligonucleotide primers used in this work are shown in 
Table S2 (Supporting Information). All PCR reactions were 
performed using total DNA of S. clavuligerus strains as 
a template in a T-gradient (Biometra, Goettingen, Germany) 
thermocycler. The PCR reaction (30 jal) was performed as 
described by Kieser and colleagues (2000), and contained 
300 ng DNA template, 0.5 mM each oligonucleotide, 28 mM 
each dGTP and dCTP, 12 mM each dATP and dTTP, 1 mM 
MgCI 2 , dimethylsulfoxide 5%, and 0.8 U Taq DNA polymer- 
ase. The amplification programme was as follows: after a 
step of 95°C for 30 s, the annealing temperature was reduced 
in a touchdown of 1 °C from 60°C to 55°C in one cycle, with an 
annealing time of 30 s; an annealing temperature of 55°C 
was used in the next 25 cycles with an extension step of 
1 min at 72°C. Quantification and purity analysis of all PCR 
products was determined using a NanoDrop ND-1000 
Spectrophotometer (Thermo Scientific, Waltham, MA, USA). 

Gene expression analysis by qRT-PCR was performed 
as previously described (Lopez-Garci'a etal., 2010). 
Streptomyces clavuligerus RNA was obtained in the same 
way as that used for microarray experiments. RNA samples 
were prepared using RNeasy mini-spin columns. The 
samples were treated with DNase I (Qiagen) and Turbo 
DNase (Ambion, Carlsbad, CA, USA) to eliminate DNA. 
Negative controls on qRT-PCR amplification (to confirm the 
absence of contaminating DNA) were carried out with each 



set of primers. The efficiency of the primers used was meas- 
ured by serial dilution of genomic DNA as template. Relative 
quantification of gene expression was performed by the AAC t 
method (Livak and Schmittgen, 2001), using the constitutive 
housekeeping gene hrdB as reference (Buttner etal., 1990). 

Luciferase assay 

For luciferase reporter analysis, promoter regions were 
amplified with primers containing Ndel and BamHI restriction 
sites (Supporting Information Table S2) to clone the promot- 
ers in the ATG codon of the luxA gene in pLUXAR-neo. 
Cultures of S. clavuligerus exconjugants harbouring the 
promoter-probe constructs were carried out in SA medium. 
Sample treatment and the luciferase assays were done as 
described by Perez-Redondo and colleagues (2012). At least 
two different cultures from the same strain were analyzed for 
luminescence production and measured in triplicate. 
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